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Abstract 

We construct explicit solutions to continuous motion of discrete plane curves described by 
a semi-discrete potential modified KdV equation. Explicit formulas in terms the r function 
are presented. Backlund transformations of the discrete curves are also discussed. We finally 
consider the continuous limit of discrete motion of discrete plane curves described by the 
discrete potential modified KdV equation to motion of smooth plane curves characterized by 
the potential modified KdV equation. 
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1 Introduction 

As is well known, many integrable partial differential equations (integrable systems) have close re- 
lationship to differential geometry. In fact, surfaces of specific curvature property in 3-dimensional 
space forms have sine-Gordon type equation as the integrability condition of surfaces. More gen- 
erally, harmonic maps of conformal 2-manifolds into semi-Riemannian symmetric spaces are con- 
structed by solutions to 2-dimensional Toda lattice equation (2DTL). 

Transformation of solutions to integrable systems have origins in classical differential geome- 
try. The Backlund transformation of the sine-Gordon equation are originally formulated as trans- 
formations of pseudo-spherical surfaces in Euclidean 3-space. 

On the other hand, substantial progress has been made in the study of discretization of inte- 
grable systems preserving "integrable structure". Motivated by extensive study on discrete inte- 
grable systems, discretizations of curves and surfaces have been recently studied actively. 
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This paper concerns with geometry of discrete curves in terms of semi-discrete integrable sys- 
tems. In [3-5], Doliwa and Santini introduced continuous motion of discrete curves in 3-sphere 
described by the Ablowitz-Ladik hierarchy [ ]. The semi-discrete potential mKdV equation was 
deduced as the simplest case. Hoffmann and Kutz [ ] introduced the notion of discrete curva- 
ture for plane discrete curves. Using the discrete curvature, they deduced the semi-discrete mKdV 
equation from continuous motion of plane discrete curves. 

In our previous works [14, 17], we have studied discrete motions of plane discrete curves in 
purely Euclidean geometric manner. The compatibility condition of a discrete motion is the dis- 
crete potential mKdV equation proposed by Hirota [ ]. In discrete differential geometric setting, 
the primal geometric object is the potential function rather than curvature (see [17]). Note that 
potential function coincides with the turning angle function in smooth curve theory. We have 
constructed explicit solutions of discrete motions of plane discrete curves in [14]. 

As a continuation of the previous works, in this paper we study continuous motions of plane 
discrete curves in terms of potential function. The purpose of the present paper is to construct 
explicit solutions to continuous motions of plane discrete curves by using the so-called r func- 
tion. Moreover we shall give Backlund transformations of continuous motions of plane discrete 
curves. The discrete curvature functions and the semi-discrete mKdV equation discussed in [13] 
are recovered from our results. 

We have been working on three categories of curves motions: (1) continuous motions of plane 
smooth curves, (2) continuous motions of plane discrete curves and (3) discrete motions of plane 
discrete curves. In this paper we investigate the relationship of these three motions, and show that 
these motions are connected by appropriate continuous limiting procedure. 

This paper is organized as follows. After recalling the requisite facts on the geometry of plane 
discrete curves and their continuous motion in Section 2, we prepare a representaion formula for 
continuous motion of plane discrete curves in terms of r function. This representation enable us to 
give explicit parametrization of motions determined by multi- soli tons as well as multi-breathers in 
the next Section 4. 

As we have mentioned before, Backlund transformation is a fundamental and effective tool 
for construction of solutions. In Section 5, we extend Backlund transformations of plane discrete 
curves studied in our previous work [14] to those of continuous motions. In particular, we give a 
new formula for Backlund transformations on the semi-discrete potential mKdV equation. 

In the final section, we shall discuss continuous limits of motions of plane discrete curves. 
More precisely, first we shall investigate continuous limits of discrete motions of plane discrete 
curves to continuous motion of those. Next we study continuous limits of continuous motions 
of plane discrete curves to continuous motions of plane smooth curves. It should be emphasized 
that these limiting procedure preserve solutions of equations. More precisely, we shall show that 
these limiting procedure preserve soliton type solutions. This is confirmed by careful analysis of r 
functions. Appendix will be devoted to detailed computations of bilinear equations for our use. 

In a separate publication [ ], we study discrete hodograph transformations and apply those 
to obtain discretizations of some integrable systems associated with continuous motions of plane 
smooth curves. 
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2 Continuous Motion of Plane Discrete Curves 



We start with the following definition. 

Z — » R 2 ; / i— > y/ is said to be a discrete curve of constant segment 



Definition 2.1 A map y 

length e if 



7m - Ji 



= 1. 



(2.1) 



We introduce the angle function ipi of a discrete curve y by 

7m ~ 7 1 



COSl/r, 

sini/f/ 



A discrete curve y satisfies 



7m ~ 7 1 



R(K,) 



7 1 ~ 7/-i 



€ 6 

for Kj = - , where R(K[) denotes the rotation matrix given by 



R(K,) 



cos Ki -sin Ki 
sin Ki cos Ki 



(2.2) 



(2.3) 



(2.4) 



We consider the following motion of discrete curves 

1 



dji _ 1 R (_ ^/\ 7/+i -7/ 
<fc cos^ V 2 



(2.5) 



Then from the compatibility condition of (2.3) and (2.5), there exists a potential function Q\ such 
that 

^/ = — j — ' ' = ? ' 

and it follows that from the isoperimetric condition (2.1) that 6 t satisfies 



d0i 
ds 



2 / - 

- tan 

6 I 4 



(2.7) 



Equation (2.7) is called the semi-discrete potential modified KdV (mKdV) equation. 

Hoffmann and Kutz [12, 13] introduced (2.5) as the edge tangential flow of discrete plane 
curves, which was deduced by discretizing the curvature function of motion of plane smooth 
curves. We note that Doliwa and Santini formulated in [3-5] the integrable motion of discrete 
curves in 3-sphere described by the Ablowitz-Ladik hierarchy, where the semi-discrete potential 
mKdV equation (2.7) arises as the simplest case. Their formulation includes the motion of plane 
curves as a limiting case. 
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3 Representation Formula in terms of r Function 

In this section, we present a representation formula for curve motions in terms of r function. We 
also give explicit r functions which correspond to soli ton and breather solutions. 

Let T[ = Ti(s;y) be a complex function dependent on the discrete variable / and two continuous 
variables s and y, satisfying the following system of bilinear equations: 



D s Ti'T* = — (r^jTz+i -T* +1 T/_i), 

T / T / = \ {tU t m +t; +1 t,_i), 



(3.1) 
(3.2) 



-D s D y TrT, = -t* 1+1 t]_ x , 
D y r l+l -ti = -eT* M T* t . 



(3.3) 
(3.4) 



Here, * denotes the complex conjugate, D s , D y are the Hirota's bilinear differential operators (D- 
operators) defined by 



D\D y f-g=(d s - d s ,)\d y - d y ) j f(s,y)g(s',y')\ s=s>y=f 



(3.5) 



We refer to [9] for calculus of D-operators. The functions satisfying the bilinear equations are 
called the t functions. 

Theorem 3.1 Let t\ be a solution to eqs.(3.1 )-(3.4). Define a real function Oi(s;y) and an K 2 - 
valued function yi(s; y) by 

2 



0,(s;y) :-- 



!og— > 



1 



- (log T(r\) y 



!og — 



2V-f 



(3.6) 



(3.7) 



Then for any s,y £R and I € Z, the functions 8i = 6i(s;y) and ji = yi(s;y) satisfy (2.1), (2.3) (2.5) 
and (2.7). 



Proof. Express ji = f (X/, Y[). From (3.4) and its complex conjugate we have 



log 

Tl 



-e 



T/+1T/' 



-6 



Tl+lTi 



(3.8) 



'/+!'/ 



Adding these two equations we obtain 

(logT M T* l+1 ) y - (log TiT]) y = -e 

which yields 



r ;+i T / T i+i T i 



Xi+i - Xi 



cos ip} 



1 

4 



log 



T/+1T/ T M T { 



w+r/ 



(3.9) 



(3.10) 



Subtracting the second equation from the first equation in eq.(3.8) we have 

Y M - Y, 



Therefore we obtain 



7m ~ 7 1 



= sini/f/. 



COSl/r, 

sin iff i 



(3.11) 



which gives eq.(2.1). Next, from eq.(3.1 1) we see that 

7m ~ 7 1 V( . , > 7i ~ 7m . . Om - di-i 

=i?((/r/-l/r/_l) , iA,-i/f,_i= = K h (3.12) 

e e 2 

which is nothing but eq.(2.3). In order to show (2.5), we identify K 2 as C. Then by using (2.2) and 
(2.6), we see that (2.5) is rewritten as 



Ki . -V^T^ ~7i yzr W^lgbi 
cos — 7i = e v 2 = e v 4 

2 6 



(3.13) 



We have 



K, 1 
cos — = — 
2 2 L 



e v i 4 + e y i 4 



T ?+l T ?-l 
T* l+l Ti-l 



1/2 



+ 



T /+1 T '-1 

T '+i T ;-i 



1/2 



1 T /+l T /-i + T /-1 T /+1 

2 



|T/+lT/_i| 



where we have used (3.2). Noticing that 

7i = x l+ V=Ty z = -(io g T;) y , 

the left hand side of (3.13) can be rewritten by using (3.3) as 

T,T* lD s D y T*-T* T M Ti-i T[ 



(3.14) 



(3.15) 



Ki dji TiT* t 

cos - -f = : —. X (-1) logr ) 



|T/+lT/_ 



(if 



|T/+lT/_i| T, 



which implies (3.13). Finally, the semi-discrete potential mKdV equation (2.7) can be derived by 
dividing (3.1) by (3.2). □ 



4 Explicit Solutions 

We now present explicit formulas for the r function which correspond to multi-soliton and multi- 
breather solutions to the bilinear equations, respectively. 

Theorem 4.1 For N e Z>o, consider the r function 

T l (s;y) = exp[-(s + el)y] del (./f, ). . | y . (4.1) 

ff = atfil - e Pi )- l e^ S+ » y + M-PiT^ + ep^'e'^^, (4.2) 
where at, /?,■ and pi (i = 1, . . . , N) are parameters. 



5 



(1) Choosing the parameters as 

p u a, 6 R, ft e V^Ir (/ = 1, . . . , N), (4.3) 
then Ti satisfies the bilinear equations (3.1 )-(3.4). This gives the N-soliton solution to (2.7). 

(2) Taking N = 2M, and choosing the parameters as 

p h a h Pi e C (/ = 1, . . . , 2M), p 2r = p\ r . x (r = 1, . . . , M), 

(4.4) 

a 2r = 0? lr _ X , for = ~P\ r -\ (? = 1, . . . , M), 

then T/ satisfies the bilinear equations (3.1)-(3.4). This gives the M-breather solution to 
(2.7). 

In order to prove Theorem 4.1, we first consider the following "generic" r function and system 
of bilinear equations. Then Theorem 4.1 is derived by applying the reduction. 

Proposition 4.2 Let cr k lm = cr k lm (u, v; y) is a function depending on three discrete independent 
variables k,l,m e Z and three continuous independent variables u, v, y e R defined by 

o^Ju^y) = det^^m)). . =j ^ (4.5) 

f®(l, m) = a iP k (l - apiY\\ - bp t y m e^i + ^ + n y +p i q k i {\ - aq t )~ l {\ - hq^e^**^*** (4.6) 

where a, b, or,-, /?,, p,- andqi (i = 1, . . . , N) are parameters. Then o~ k lm satisfies the following bilinear 
equations: 

(A, - 1) a£? • o-l = -o- k l+hm o- k :l m , (4.7) 
(D v - 1) o-U ■ o-l = -< + i<-i> (4.8) 
bo&Ki* ~ "<L< + i + (fl ~ Mtl«,Am = 0, (4.9) 
l -D u D y < • < = a(<) 2 - art\y-{ m , (4.10) 

^D v D y < • Cr k m = b(cr k n f - bcr^O-lU, (4.11) 

(D y - a) oj +hm ■ < = -aoH^oft. (4.12) 

Proof of Theorem 4.1 We show that Theorem 4.1 holds from Proposition 4.2. We impose the 
reduction conditions on o~ k as 

l,m 

^Vi, m+ i = B^ m , (4.13) 
< m l =Co-* k , CgR (4.14) 
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where B, C are constants. Then putting b = -a, the bilinear equations (4.7)-(4.12) are reduced to 

(D u - 1) &\ ■ o-i = -o- M o-\_ x , (4.15) 

(A- - 1) cr* ■ o-/ = -cr,^cr* +l , (4.16) 

0"/_iCn+i + cr^ +1 cr M - 2cr*cr/ = 0, (4.17) 
1 2 

-D u D y cr r ai = a(cr,) - aa* M cr\_ x , (4.18) 
1 

-D v D y a r o-i = -a(cri) + acr^cr]^, (4.19) 

(p y - a) cr /+1 • cr, = -acr* M cr), (4.20) 

respectively. Here we have used (4.13) and (4.14) to eliminate the m- and ^-dependence, respec- 
tively, and denoted cr k lm = cr,. We next consider the specialization of continuous independent 
variables 

u = cs, v = -cs, cel. (4.21) 
Then, subtracting (4.16) from (4.15) we have 

D s cr*, -cri = c (cr^icr* l+l - (Ti+iCr*^) . (4.22) 

Similarly, we get from (4.18) and (4.19) 

D s D y cr r cr, = 4ac [(cr,) 2 - crUcr* M ) . (4.23) 

Putting 

a = 6, c = (4.24) 

and introducing t/ by 

Tl = e~ (s+d) >cr h (4.25) 

the bilinear equations (4.22), (4.17), (4.23), (4.20) are reduced to (3.1), (3.2), (3.3), (3.4), re- 
spectively. Let us next realize the reduction conditions (4.13) and (4.14) by imposing suitable 
restriction on parameters of solution. We put 

qt = -Pi (i=l,...,N), b = -a. (4.26) 

Then it is easy to verify that the entries of the determinant satisfy 

/f (Z + l,m + 1) = — ^— r yf (Z,m), (4.27) 
1 - a z pj 

so that the condition (4.13) is realized as 

<w = fl r^v ^ (4 - 28) 

;=l w 
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As for the condition (4.14), we have to consider the cases (1) and (2) in Theorem 4.1 separately: 
Case (1). We impose the condition (4.3). Then we see that 



f® l (l,m) = p i f«>*(l,m), 



■CO*/ 



and so 



a k ^=Cal nv C = Y\ Pi eR. 

i=l 

Case (2). We impose the condition (4.4). Then we see that 

/£?(/, m) = p* 2r _jl lr - l) \l, m), jf:; l \l, m) = p\ r ff r) \l m), 

and so 

M 



°%=Co%, C=(-l) M Y\\P2r\ 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



r=\ 



Finally, putting m = without loss of generality and applying the specialization (4.21) and 
(4.24), (4.6) is rewritten as 

/® (Z,0) = a iP k {\ - epiT'e^^^y^' +fr(- Pi ) k (l + epiT'e^^"^^^ 



atPid - epd l e 



+ j3 i (-p i ) K a + ep i y l e 



which is equivalent to (4.2). Therefore we have derived Theorem 4.1 from Proposition 4.2. □ 

The bilinear equations in Proposition 4.2 are reduced to the quadratic identities of determinants 
(Pliicker relations). In particular, (4.9) and (4.12) have already appeared in [14]. Moreover, by the 
symmetry between the set of variables (/, u) and (m, v) in a k lm , it suffices to show only (4.7) and 
(4.10). These bilinear equations will be proved in the Appendix. 

Remark 4.3 In the r function in Theorem 4.1, the parameter dependence of the time evolution 
in entries of the Casorati determinant have singularities different from and oo. These types of 
singularities can be seen in the solutions of equation of principal chiral fields, i.e., harmonic maps 
of conformal 2-manifolds into compact Lie groups [15,24,25] and Maxwell-Bloch equation [16]. 

Remark 4.4 By introducing u\ := the semi-discrete potential mKdV equation (2.7) can be 
transformed to the semi-discrete mKdV equation 



du 



— = [I +U 1 J(U M -M/_i), 



(4.33) 



where we put s = 2es' for convenience. An auxiliary linear problem for (4.33) is given by [3] 



= 



1 



1 + uf 



A A x u\ 
-Aw A~ l 



A 1 -A 



ui + A ui-i 



—U[ — A U[-\ 



A L -A- 



O,. (4.34) 
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Apparently, the dispersion relation suggested from the linear problem is different from the one in 
Theorem 4.1. However, putting 

1 X] - 1 

Pi = €aJTT (435) 

in (4.2), then can be rewritten as 

in which the dispersion relation with respect to I and s' is consistent with (4.34). We have chosen 
the parametrization as in (4.2) so that the continuous limits explained in Section 6 become simpler. 



5 Backlund Transformations 

In this section we discuss the Backlund transformation of the continuous motion of plane discrete 
curves. The Backlund transformation of the plane discrete curves has already been formulated 

in [14]: 

Proposition 5.1 Let ji be a discrete curve of segment length e. Let 0[ be the potential function 
defined by 



7/+1 - 7n 



COS)/// 

sin </// 



0i+ i + Oi 



(5.1) 



For a nonzero constant A, take a solution 6 n to the following equation 



tan 



(0 



7+1 - m 



tan 



1~ 6 



then 



~ 1 D ( 01 ~ 0M 



01 - 0l+i 



7i+\ ~ 7i 



(5.2) 



(5.3) 



is a discrete curve with the potential function 6\. 

We next extend the Backlund transformation to that of motion of discrete curves. In order to do 
so, we first present the Backlund transformation to the semi-discrete potential mKdV equation: 

Lemma 5.2 Let 0i be a solution to the semi-discrete potential mKdV equation (2.7). A function 9[ 
satisfying the following system of equations 



e tan = — h e I tan 

U 4 \A 



4cos : 



01-0M \ A 



4 cos 2 







7+1 - 01 



01 - 01+1 01+i - 01 

tan ■ + tan ■ , 



4 4 

gives another solution to eq.(2.7). We call 0\ a Backlund transform of0\. 



(5.4) 
(5.5) 
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Proof. First compute addition of (5.5)/_i and the derivative of (5.4)/_i . Then, by using (5.4), 
eliminate A from this equation and (5.5) respectively. Adding those two equations yields 



e 77/ di+i - 

- 9, cos 

2 ' 4 



9i-i . 9m 
sin ■ 



= | - ^ COS ^ 



sin ■ 



4 



9i+i - 9, 



sin 



9i+y + Oi-i - 29 j 



(5.6) 



i-i 



sin 



9i+\ + 9 1-\ - 29 1 



which implies Lemma 5.2. □ 

Proposition 5.3 Let ji be a motion of discrete curve. Take a Bdcklund transform 9i of9\ defined in 
Lemma 5.2. Then 



7i 



7l + - R 

H A 



/ 9/ - 9i+i 



7m ~ 7 1 



(5.7) 



is a motion of discrete curve with potential function 9/. We call ji a Bdcklund transform ofji. 

Proof. It suffices to show that ji satisfies eqs.(2.1), (2.3) and (2.5) with potential function 9 h but 
eqs.(2.1) and (2.3) follow from Proposition 5.1 immediately. Because the system (5.4)-(5.5) yields 



~^ + ~ 2~ ■ 



we identify R 2 with C, so that the motion ji satisfies 



- ^ 



7i=e 



7/+i - 7i ' 



l- V=T 



tan 



9m — 9, 



M 



which implies (2.5) with 2K\ = 9m - 9m- 



□ 



Remark 5.4 In [12, 13], the Backlund transformation of the motions of discrete plane curves de- 
scribed in this paper is characterized by the cross ratio of the four points ju 7m, 7t and jm being 
constant. In fact, we can verify by direct computation that for the Backlund transformation given 
in Proposition 5.3, the cross ratio of those four points is -A 2 e 2 . In case of continuous motions of 
smooth plane curves, Darboux matrices of the Backlund transformations are obtained in [ ] . Sym 
gave explicit formulas for smooth plane curves derived from multi-soliton solutions to the mKdV 
equation via the iterated Backlund transformations. 



6 Continuous Limits 

In [14], the discrete motion of discrete plane curves and the continuous motion of smooth plane 
curves have been formulated, together with the Backlund transformations and the explicit formulas 
in terms of the r functions. They are described by the discrete potential modified KdV equation and 
the potential modified KdV equation, respectively. In this section, we present the two continuous 
limits: one from the discrete motion of discrete plane curves to their continuous motion discussed 
in the preceding sections, another one from the continuous motion of discrete plane curves to the 
continuous motion of smooth plane curves. 
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We first summarize the formulations of three kinds of curve motions and explicit solutions. For 
convenience, we identify Euclidean plane K 2 with complex plane C. 

(1) Discrete motion of discrete plane curves. 

Motion of curves: 

m _ m 
I n+l in 



= h 



■V" — V" , V — 'V m 

[n+l In _ -fAK™ >n I n-l 

Cl n &n- 1 

V" +1 - V* /— y m - y m 

Y_n Y_n_ _ ^yZiw" /n+l 'n 

b m ci n 



(6.1) 
(6.2) 
(6.3) 



Here, n, m e Z denote the discrete independent variables corresponding to space and time, respec- 
tively. Moreover, a n , b m are real arbitrary functions of the indicated variables, which correspond to 
the segment length of the curves and time interval, respectively. 



Potential function: 



Compatibility condition: 



cm _ cm nm+l _ am 

U n+l U n-l w m u n °n+l 



tan 



b m ci n 



cm+ i cm \ i , I nm+ i cm 

Vn+l ~ tf n\_V m +a n (tf n ~ ff n+1 

tan 



Explicit formula in terms of r function: 



(6.4) 



(6.5) 



2 T m 

er = log — y m 



Soliton type solutions: 



t? = exp 



r 

log^ 

c j-*m 



. 2 V 3 ! 



y J 



' n-l m—l 



J]a n , + J]bJy dot . , v . 

V n' m' ) 



(6.6) 



(6.7) 



n— 1 m—1 n-l m—1 

jf = a iP ) Y\(\-a n , P fr x Yld-b^pd^e^+M-Pif Y\(\+a n , Pi y l ^(l+^p,)" 1 ^^. (6.8) 



m! 



(2) Continuous motion of discrete plane curves. 

Motion of curves: 

7/+i - Ji 



= 1, 



77+i ~ 7< _ e V^Tsr, 77 ~ 7/-i 



dji e _v_1 T y M - 



7i 



ds cos f 



(6.9) 
(6.10) 

(6.11) 
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Potential function: 



Compatibility condition: 



K,= 



dOi 2 ( 6i + i - 61-1 



- tan I 

ds e \ 

Explicit formula in terms of t function 



(6.12) 
(6.13) 



2 i Tl 
#/ = ^=log— , yi = 



-- (log TiT t ) y 



2V-F 



Soliton type solutions: 



r, = exp[-(s + e/)v] del (./;•', ) ; / | y . 



/f = a i p k i (\-ep i r l e^ Pi + /3,<- Pi ) k (l + e Pi T l i '"^ 



s- — y 



V p, 



(6.14) 

(6.15) 
(6.16) 



(3) Continuous motion of smooth plane curves. 

Motion of curves: 

\Y\ = i, 



— y = yPi K y, 

ox 



(6.17) 
(6.18) 

(6.19) 



Here y = y{x, t) e R 2 C is arc-length parametrized curve, x and t denote arc-length and time, 
respectively, and ' = d x . Moreover, k = k(x, t) is the curvature. 



Potential function: 
Compatibility condition: 



1 



(6.20) 
(6.21) 



Explicit formula in terms of r function: 



2 , T 

6, = -f= 1o §— ' y = 



-- (log TT*) y 
Pi V T* />• J 



2^ 



Soliton type solutions: 



r = <T^ det(/ (!) ,) 



(6.22) 

(6.23) 
(6.24) 
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Theorem 6.1 



(1) Putting 

a + b a-b 

a n = a (const.), b m = b (const), 5 = — - — , e = — - — 



s , 

- = n + m, I = n — m, 

6 



2 - (6.25) 



and taking the limit 8 — > 0, the discrete motion of discrete plane curves yields the continuous 
motion of discrete plane curves. 



(2) Putting 



x = el + s, t = - — s, (6.26) 
6 

and taking the limit e — > 0, the continuous motion of discrete plane curves yields the contin- 
uous motion of smooth plane curves. 

Theorem 6.1 can be verified by tedious but straightforward calculations. In fact, the statement 
(1) can be checked by substituting the parametrization (6.25) into (6. 1)— (6.8), expanding in terms 
of powers of 6 and taking the limit £ — > 0. The statement (2) is also checked by a similar manner. 
We note that the limiting procedures presented in (6.25) and (6.26) have been obtained in [7] 
and [ ] on the level of the equations for 6. Theorem 6.1 claims that the procedure applies to the 
curve motions and solutions. Also, it should be noted that limiting procedure also applies to the 
Backlund transformations. 

In order to demonstrate the calculation, we here discuss the limits of the r functions corre- 
sponding to the soliton type solutions. Substituting (6.25) into (6.8), we have 



(1 _ ap . r n (l _ bpi ym = (1 _ aA .)-Kf+0(l - %)4(!-0 = c -*Jog[l-*W^»tf] ( I 6Pi 5 ll 

\ \ + epi-Spt, 

Noticing that 

log [ 1 - 15 Pi + (6 2 - e 2 )p]} = log coi -^6 + 0(6 2 ), oj t = 1 - e 1 p], 

<X>i 

we get 



(1 - a Pi y n (\ - b Pi T m ~ e~— s x e~< s ' 1 



1 + epi 



Similarly, we have 



log", Pi . / 1 + 6 Di 

(1 + a Pi y n (l + b Pi )~ m ~ e~— 1 x e~~ s 1 1 



1 - epi 



Therefore /, (!) yields 



e~—\\ - e 2 p]Y 



a iP 1 (1 - e Pi y e^ y +#(-#)*(! + e Pi T l <r^ y 



(6.27) 
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as 8 ~ 0. The prefactors of the entries in (6.27) can be factored out of the determinant, and it is 
easily seen that the overall factor does not affect the solutions, namely, if we remove overall factor 
from the r functions, it gives the same 6 and y, as seen from (6.6). This implies that the determinant 
in (6.7) yields that in (6.15) up to this trivial multiplicative factor. Also, the exponential factor in 
(6.6) becomes that in (6.14) under the parametrization (6.25). Therefore, we have shown that (6.7) 
is reduced to (6.15) as 8 — > 0. 

Similarly, substituting (6.26) into (6.16), we have 



(1 - epd l e'- e2p r = exp 



and 



- + ^|log(l -epi) + 



l-6 2 p 



2„2 



= exp 



-^-t + ( Pi x - 4p]t) + 0(e) 



(l + e Pi T l e l -' 
from which we obtain as e ~ 



< <.■ 



exp 



-^-t-(p i x-4p 3 i t) + 0(e) 



fi 



(0 



3/7 



k PiX-4ph+^-y , D / \k -PiX+4-ph-^-y 



(6.28) 



The prefactor in (6.28) does not affect the solutions. Also, the exponential factor in (6.14) becomes 
that in (6.22) under the parametrization (6.26). Therefore, we have shown that (6.15) is reduced to 
(6.23) as e -» 0. 

A Derivation of bilinear equations (4.7) and (4.10) 

In this appendix we prove Proposition 4.2. As mentioned in Section 3, it suffices to show that the 
r function given in (4.5) and (4.6) actually satisfies the bilinear equations (4.7) and (4.10). 



A.l Equation (4.7) 

We define the r function cr k lm {u, v; y) by 

(rj>,v;j) = detf/^am)). . =J n , 

where the entries of determinant satisfy the linear relations 

fl\l, m) - /®(/ - 1, m) /»(/, m) - /f (/, m - 1) 



a ' b 

»•>// ,„> _ , i a AO st ^ _ ADn _ , i% * AOs, _^ _ AO 



(A.1) 



(A.2) 



duf?>(!,m) = f?(l + l,m), d v f^(l,m) = f*>(l,m+ 1), d y f^(l,m) = /^(/,m). (A.3) 



Note that f^\l, rri) given in (4.6) satisfy the above relations. In order to prove (4.7), it is convenient 
to consider pf defined by 



p\ >, v; y) = det (jf (Z - j + 1 , m)) . , = ^ 



(A.4) 
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instead of <x* . Here, of and pf are related as 

l.m ' Lm ' Lrn 



n k 

Pirn 



(-a) 



JV(JV-l)/2 _£ 



/,m' 



(A.5) 



which can be easily verified by manipulating the columns of determinant with the first equation in 
(A. 2). We also introduce a notation 



tf i k ... N — 2 k N — V 



fl l \l-j,m) 
jf\l-j,m) 



fl N \l-j,m) 



(A.6) 



It is possible to reduce (4.7) to one of the Pliicker relations which are quadratic quadratic identities 
of determinants whose columns are appropriately shifted. To this end, we construct such formulas 
that express the determinants in the Pliicker relations in terms of derivative or shift of discrete 
variable of p k lm {u, v;y) by using the linear relations of the entries. For details of the technique, we 
refer to [9,18-21]. 



Lemma A.l The following formulas hold: 



d up\ m = 1 "1,1, 


■ ,N- 


-2,N- 


-H, 


plt=a N ~ l |0,1,- 


■■ ,N 


-2,N 


-l k ~ l \, 


ptn=a N ~ l |0,1,- 


■■ ,N 


-2,N 


-2 k ~ l \, 


a - 1) pj; 1 = a N ' 


i | _ 


1,1, 


■ ,N-2,N-l k 



k-l\ 



(A.l) 
(A.8) 
(A.9) 
(A. 10) 



Note that the superscripts of column vectors are shown only when k is shifted for notational sim- 
plicity. 

Proof. Equation (A. 7) follows from the differential rule of determinants and the fist equation of 
(A.2). Next, applying the first equation of the difference rule (A. 3) to the first column of p)z}, we 
have 



• ,N-\ k ~ l \ = I 0* _1 - I* -1 ,!* -1 ,-- - .iV-l* -1 ! = a I 0*,1* -1 , 



Pirn = |° , , , , 

Repeating this procedure for the j'-th column (j = 2, 3, . . . ,N - 1), we get 

Pu = a 



,N-1 



k-\\ 



- v |0*, ,N-2 k ,N-l k ~ l \ 



which is (A.8). Applying (A.2) to the N-th column of (A.8), we obtain (A.9). 
Finally, differentiating (A.9) by u yields 



°u Plm ~ a 



JV-1 



-l k ,l k , ••• ,N -2 k ,N -1^ +a N ~ l |0*,1*,--- ,N -2 k ,N -2 k ~ l \ 
-l k ,l k ,--- ,N-2 k ,N-l k - 1 \+p k ;^ 
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ft- II 



which is equivalent to (A. 10). Thus we have proved Lemma A.l. □ 

Now consider the Pliicker relation (see, for example, [21]) 

= | - 1, 0, 1, • • • , N - 2\ x 1 1, • • • , N - 2, N - 1, N - 1 
+ | 0, 1, - - - ,N-2,N-l\x\ -1,1, ,N-2,N-l k ~ l \ 

— | o, i, • - • ,iy-2,iy-i 4_1 | x i ,iv-2,iv-i|. 

(A.l 1) is rewritten by using Lemma A.l as 



(A. 11) 



= Plun x a 



-(N-l)k-l 



-(AT-1), 



PT-im + Pl m x a— (d u - 1) pf" 1 - a-^ptt x d u pi, 



(A.12) 



= fl-W-« [(D u - 1) pf" 1 -pj w +p k Mtn fit{ m ] , 
which implies (4.7). 

A.2 Equation (4.10) 

We derive (4.10) from (4.7) and (4.12). We first introduce F k m by the subtraction of the right hand 
side of (4.10) from the left hand side 



(A.13) 



and consider 



p ._ fk t^k-lS} _ pk-l / i \ 2 
' ■- r l,m\yijn ) r l,m \yi,m) 



«) 2 



Equation (A. 14) can be rewritten as 



\DuD y o 3 ^ • a k ^ - a(ojJf + acr k +hm cr k :l m 



(A. 14) 



P = D y (D x rf m ■ cr k J) ■ a-\y^ + ao*:l m oj:l m o*^(j*^ - aoj^cr^oj^, (A.15) 



where we have used the exchange formula of the D-operator [9] 

(AA f-f)g 2 -f (D u D y g-g) = 2D y (D u f ■ g) ■ fg, 



(A.16) 



for arbitrary functions / and g. We manipulate the first term of (A.15) as follows. Using (4.7) and 
noticing D y f • f = 0, we have 

„ / „ £ k—l\ k k—\ r\ f k—lk k k—1 \ k k—\ t\ k k—\ k k—l 

D y [ D x Vim ■ °l,m ) ■ ^InFUm = D > \-°~lm ^l.m + °t + l,m°/-l, m J ' ^nPljn = A> °l + \,rrPl-\jn ' rfmrfm ■ 

Then applying another exchange formula 

D y a/3- yd = (p y a-y)/36 + (p y yS • s) ay, (A. 1 7) 

16 



for arbitrary functions a, /?, y, 5, we get 

Dy cr l+l,m oJ l-ljn ' °^m cr tii = fay °~l+\,m ' °~l,m) <J f-l,ro £r tn + fay ^f-l,™ ' °"tn ) fJ f+l,m Cr t»j 

where we have used (4.12) in the second equality. Substituting the above result into (A. 15), we see 
that P = 0. Therefore, it follows from (A. 14) that 

l -D u D y a\ m • erf m - a(o-[J + a</, m ^ m = A(u, y, /)(<) 2 , (A. 1 8) 

where A(u,y, I) is an arbitrary function. Since a k lm = 1 (the case of Af = 0) satisfies (4.7) and 
(4.12), it should satisfy (A. 18) as well. Therefore we see that A must be 0, which implies (4.10). 
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